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SECTION  I 

PROGRAM  DESCRIPTION 


The  program  "OPTIM"  Is  an  optimization  program  designed  to  run  on  a 
small  microcomputer.  The  listing  presented  here  Is  for  a  Sharp  PC-1500 
with  an  8K  memory  expansion.  The  version  of  Basic  Is  standard  enough 
that  translations  to  other  microcomputers  Is  possible. 

The  program  will  find  a  minimum  thickness  lamlmate  which  will  not 
fall  under  any  of  the  load  conditions  entered.  Ply  orientations  are 
chosen  by  the  user.  The  program's  capability  In  handling  multiple. 
Independent,  loads  could  be  useful  for  loads  which  change  with  time  or 
for  situations  where  there  Is  uncertainty  In  calculating  the  loads.  As 
the  program  Is  currently  dimensioned,  four  Independent  load  combinations 
and  6-ply  orientations  can  be  entered. 

Only  point  stresses  are  considered,  thus  the  program  optimizes  the 
laminate  only  at  one  point  In  the  structure.  Furthermore,  the  program 
assumes  In-plane  loads  only  and  no  out  of  plane  deflections.  This 
implies  a  symmetric  laminate,  but  stacking  sequence  Is  not  a  factor  In 
the  program.  The  layer  thicknesses  generated  by  the  program  are  the 
total  and  must  be  divided  by  2  to  get  the  halves  of  a  symmetrical 
laminate. 

No  knowledge  of  optimization  techniques  Is  needed  to  run  the  program 
and  very  little  knowledge  of  laminate  plate  theory  Is  needed.  In 
addition,  material  properltles  for  three  common  advanced  composites  are 
stored  In  the  program  and  new  composite  materials  can  be  added  by 
Inputting  their  material  constants  Into  the  program  directly. 
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SECTION  II 

GENERAL  INSTRUCTIONS 


The  Sharp  PC-1500  Manual  Includes  tape  loading  Instructions.  Because 
"OPTIN'1  takes  so  long  to  read  (approximately  5  minutes).  It's  a  good  Idea 
to  test  tape  player  volume  level  with  a  short  one  or  two  line  program  to 
see  If  all  Is  well.  Load  the  tape  using  "OPTIN"  as  the  name.  An  example 
of  the  display  end  appropriate  responses  are  given  In  this  report. 

The  program  Is  started  In  RUN  mode  by  Instruction  RUN  (press  the  keys 
[R]  [U]  [N]  and  [ENTER]).  User  Is  guided  through  the  program  by  simple 
questions.  The  user  types  the  chosen  answer  and  presses  the  [ENTER] 
key.  At  the  completion  of  the  routine,  after  all  results  have  been 
displayed,  the  program  will  restart  Itself. 

Run  times  can  be  quite  long.  They  range  from  a  minute  for  a  2-layer 
laminate,  to  10  minutes  for  a  6-layer  laminate  subject  to  multiple  loads. 

The  material  constants  are  stored  In  the  program.  When  "MATERIAL 
NO.*"  Is  asked.  Input  the  material  number; 


1.. ..T300/5208 

2.. .. B0R0N/5508 

3.. .. AS/3501 


To  add  other  materials  (up  to  6  more),  user  should  type  (In  PRO  mode) 
for  a  material  with  the  material  number  N, 


16N?:  DATA  "name  of  material",  E  ,  E  ,  Vy,  E  ,  h  ,  X,  X',  Y, 

A  A  5  U 

Y',  S  :  RETURN. 


An  example  of  adding  a  new  material  (Aluminum.... Material  number  4)  Is  as 
follows: 


Set  PRO  Mode  and  type, 

1640  DATA  "ALUMINIUM",  69,  69,  .3,  26.538,  125E-6,  400  ,  400,  400,  230 
RETURN 

and  press  ENTER  key. 

The  unit  of  engineering  constants  are  In  GPa  and  thickness  of  single  ply 


hQ  Is  In  meters. 


SECTION  III 


METHOD 


The  goal  Is  to  minimize  the  total  thickness  of  a  composite  laminate 
subject  to  failure  constraints  under  static  loads.  Specifically 


k  h,  =min.  where  L  =  number  of  layers 

k=l  k 


subject  to  h  >  0  and  G  Gj9kMN)e(N)+  Q(Qk^e(N)  -  i  <  o  where  hk  Is 
the  total  thickness  of  all  the  plies  at  the  k'th  orientation  (which  will 
be  referred  to  as  a  "layer"  In  this  report).  The  failure  criteria  Is  a 
first  ply  failure  based  on  the  Tsal-Mu  tensor  criteria  In  strain  space. 
The  G's  are  transformed  to  the  laminate  axis  from  the  k'th  layers 
orientation.  The  strains  are  associated  with  the  N'th  loading 
combination.  This  distinction  Is  made  since  more  than  one  Independent 
loading  may  be  considered.  For  the  definition  of  the  G's  In  terms  of 
experimental  strength  data,  see  reference  1. 

Stacking  sequence  Is  not  Included  In  this  formulation,  and  the 
laminate  Is  assumed  not  to  bend  or  warp.  Therefore,  strains  and  loads 
are  related  by 


!$  =  |A|  t 

The  optimization  method  applied  Is  a  modification  of  the  method  of 
feasible  directions.  The  method  can  be  demonstrated  graphically  with 
2-dimensions,  l.e.  two  layers.  In  figure  1  the  two  equalities 

Gij0)vj  *  g!°S  - 1  ■  ° 


G<?°>eiE. 
U  1 7  J 


♦  G<9°>  e. 
1  1 


1  =  0 
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have  been  plotted  as  functions  of  h^  and  h^®0^  for  the  single 
loading  condition  shown.  Any  point  above  and  to  the  right  of  these  two 
curves  Is  feasible,  that  Is,  failure  will  not  occur.  Points  to  the  left 
and  below  the  curves  are  Infeasible.  Because  our  objective  function  (the 
sum  of  the  layer  thicknesses)  Is  linear,  the  optimum  point  will  He  on 
one  of  these  curves  or  the  Intersection  of  multiple  curves. 

The  program  starts  by  finding  an  Initial  feasible  point  (A)  which 
lies  on  a  constraint  curve  farthest  from  the  origin  on  the  line 
hCSO]*^].  ^  distance  from  the  origin  Is  calculated  using  a 
strain  ratio  method.  Along  any  vector  which  passes  through  the  origin 


s/s 


o 


where  S  Is  a  scalar  distance,  and 


s„  =  [  Uw/2 


along  this  vector,  strain  can  be  found  using 


where  e^is  a  component  of  laminate  strain  evaluated  at  SQ. 
Substituting  Into  the  failure  criteria  we  have 


To  ensure  the  calculated  point  lies  sllghty  In  the  feasible  region 
despite  any  numerical  error,  the  program  sets  this  function  equal  to  the 
negative  of  a  small  number  lE(-6)  rather  than  zero.  Solving  this 

equation  for  positive  S  we  have 

-B  +  /(B2-4AC) 


where 


<rk'£?s0 


A=1-lE-6 


B"i=l 


3  3  (q  )  °  0  o 

C=iil  jSj  -gWc-e^ 

If  SQ  lies  in  the  feasible  region  we  solve  the  above  equation  for  each 
layer  and  each  load  combination  then  take  the  smallest  resulting  S  as  the 
one  that  defines  the  boundary  of  the  feasible  region. 

The  next  step  in  the  optimization  procedure  Is  to  establish  a 
direction  vector  which  will  point  away  from  the  constraint  A  lies  on  and 
Is  parallel  to  the  plane  defined  by  £  h^constant.  In  figure  1,  this 
direction  Is  shown  as  Z.  Finding  Z  first  requires  calculation  of  the 
gradient  of  the  active  constraint  evaluated  at  A.  Let 


C  =  G<ek>e(N>r  <N>  +  G(ekMN) 
Ck,N  bijk  ei  j  i 


1 


where  k  and  N  correspond  to  the  layer  and  load  combination  of  the  active 

\ 

constraint.  A  constraint  is  considered  active  if  > -E2  where 
E2"lE-5  is  a  small  number.  Note  that  more  than  one  constraint  may  be 
active.  The  gradient  Is  then  given  by 


$  C, 


k,N 


*  E  [G^k* 
L*1 


3ei 


(N) 


.(N)  ■  AH) 


3h, 


+  £i 


_ J. 


(N) 


ah, 


)  + 


9ei 


(N) 


3h, 


where  hL  Is  a  unit  vector, 
with  the  basic  equation 


To  find  the  parti als  of  strain  we  start 

n-IaI  * 


O.-i-  |A|  e  +  A  -2-e 


iL  =  _| A-1 1  ~  |A|  e 


and 

Re,-> 

«n 

(0|) 

*12 

(e.) 

*13 

—  i*i  - 

3hi 

(0,) 

(V 

*22 

(e.) 

Q23 

(0,) 

«13 

(<v 

Q23 

(0.0 

Q33 

The  gradient  vector  Is  normalized  to  unit  length.  If  more  than  one 
constraint  Is  active,  the  normalized  gradients  are  summed  together  and 
the  sum  Is  then  normalized  to  one.  The  negative  of  the  gradient  will 
point  away  from  the  constraint.  Into  the  feasible  region.  This  vector  Is 

A 

now  projected  onto  the  plane  defined  by  the  unit  normal  n,  where 


A  1  L  A 


The  projection  can  be  made  with  a  double  cross  product 

A  A 

Z  *  n  x  (-Vc  x  n)  . 

With  a  vector  Identity,  this  can  be  rewritten  as 

_  A  A  - 

Z  *  (Vc  •  n)n  -  Vc  . 

-*■ 

Finally,  Z  Is  also  normalized  to  unit  length. 

Along  Z,  another  constraint  will  eventually  be  reached  (point  B  in 
Figure  1).  The  point  Is  found  Iteratively  by  a  bisection  technique. 

Since  the  bisection  method  Is  very  time  consuming,  the  constraint  line  Is 
only  found  within  a  relatively  large  error  band.  What  we  are  really 
Interested  In  Is  a  point  approximately  midway  between  A  and  B,  which  Is  C 
In  the  figure.  From  point  C,  the  strain  ratio  technique  Is  used  to 
analytically  calculate  D.  Starting  at  D,  the  entire  procedure  repeats. 
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The  program  terminates  when  the  distance  AS  or  CD  Is  small  (say  1/10  of  a 
ply  thickness)  or  the  magnitude  of  Z  before  normalization  Is  very  small 

A 

(Implying  n  and  vc  are  almost  parallel). 

In  some  cases,  h^O  constraint  may  be  reached.  When  this  happens, 
that  orientation  Is  completely  dropped  from  further  calculations.  Thus, 
the  constraints  associated  with  a  zero  thickness  layer  cannot  affect  the 
results.  Once  an  orientation  reaches  zero  thickness.  It  Is  never 
reinstated  In  later  Iterations. 

Figure  1  shows  a  case  where  the  program  reaches  the  Intersection  of 
two  constraints.  However,  simultaneous  failure  should  not  be  considered 
a  criteria  for  optimization.  The  constraint  line  for  the  +45°  layer  Is 
completely  In  the  Infeasible  region.  The  line  h^-45^  +  h^-"45^  *  const, 
has  been  Included  to  show  that  point  D  Is  the  minimum  thickness  (See 
Figure  2). 
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Figure  1.  Constraint  Surfaces  and  Optimization 

Trajectory  for  0/90  under  Biaxial  Load. 


Figure  2.  Constraint  Surfaces  and  Optimization  Trajectory 
for  ±  45  Laminate  Under  Pure  Shear. 


APPENDIX  B 


KEY  OPERATION  PROCEDURE 


-  display _ 

INPUT 

PRINT  OUT  *  REMARKS 

MATERIAL  NUMBER  - 

RUN  [ENTER] 

1  [ENTER] 

T300/5208 

NUMBER  OF  PLY  ANGLES  ? 

PLY  ANGLE-1  ? 

2  [ENTER] 

0  [ENTER] 

NO. OF  PLY  ANGLE-  2 

ANGLE-  0.0 

PLY  ANGLE-2  ? 

90  [ENTER] 

ANGLE=  90.0 

NUMBER  OF  INDEPENDENT 

LOADINGS 

2  [ENTER] 

LOADING-1 

Nx=? 

4  [ENTER] 

NO. OF  LOADINGS=  2 

LOADING  UNIT  is 

N2»? 

1  [ENTER] 

in  riN/m 

Ng-? 

0  [ENTER] 

LOADING-1 

Nl=  4.000 

N2=  1.000 

LOADING-2 

N6=  0.000 

Ni»? 

1.25  [ENTER] 

LOADING-2 

n2»? 

3.25  [ENTER] 

Nl=  1.250 

N2=  3.250 

N6=  i . 299 

n6«? 

1.299  [ENTER] 

I V.  -Xt- 


PRINT  OUT  l  REMARKS 


ANGLE=  0.0 
RAT  J  0=  0.304 

#  PL1ES=  55.63 

ANGLE=  90.0’ 

AT  10=  0.695 

PL  I ES=  126.94 

STRENGTH  RATIOS 

.1  =ULT IriATE  STRAIN 
AND  > ] =SAFE 

LOADING-1 
PLY  RATIO 

>3.0  2.301 

30.0  1.371 

LOADING-2 
PLY  RATIO 

0.0  I . 000 

90.0  .1.023 


—BREAK— 


APPENDIX  C 


MEMORY  CONTENTS 


MEMORY 

DESCRIPTION 

MEMORY 

DESCRIPTION 

A(3,3) 

Value  of  Thickness 

Z(6) 

Direction  Vector 

B(6,9) 

6jj  (Strength  Parameters  In 

J  Strain  Space) 

A$(6) 

Subscript  Indicators 

C(6,9) 

(Modulus  Components) 

U1—U5 

Modulus  invariants 

0(3,3) 

l  Zk 

VI— V7 

Strength  invariants 
(Strain  Space) 

E(4,3) 

Strains 

C2 

COS  20 

6(3,3) 

Strength  Parameter  Matrix 

C4 

COS  4e 

N(4,3) 

Loads 

S2 

sin  20 

P(3,3) 

Inverse  of  A  Matrix 

S4 

sin  40 

Q(3,3) 

Modulus  Matrix 

S 

Final  scalar  distance  to 
be  moved 

C$(10,2) 

List  of  active  constraints 

S1-S2 

Points  in  feasible  region 

H(6) 

Thickness  for  each  ply 

SR 

Distance  to  origin 

R(3) 

Intermediate  results 

SM 

Distance  to  first  h=Q 
constraint 

S(3) 

Strength  parameter  components 

NP 

Number  of  ply  orientations 

T(6) 

Angle  of  each  ply 

NL 

Number  of  independent 
loadings 

MEMORY  CONTENTS  (CONTINUED) 


MEMORY 

DESCRIPTION 

MEMORY 

DESCRIPTION 

W(6) 

Normalized  gradient  of  a 
constrai nt 

NC 

Number  of  active  constraints 

X(6) 

Normalized  sum  of  gradients 

IT 

Iteration  counter 

Y(3) 

Intermediate  results 

TP 

Individual  thickness 

IM 

Maximum  number  of  iterations 

_ 

GRtGS 

GP,GQ 

Strength  parameter 
in 

strain  space 

NN 

M 

1 

<. 

X 

1 

GX.GY 

mm 

Intermediate  calculations 

F$ 

Flag  to  halt  program  when  "F" 

. 

G$ 

w!n?atpH  if  a  constraint  is 

I,  J,  K, 

L 

Loop  counters 

M$ 

Material  Numbers 

P.  II 

Ply  orientation  pointers 

H$ 

Material  Names 

N 

Load  pointers 

x$ 

Subscript  indicators 

APPENDIX  D 


(0)  SAMPLE  PROBLEMS 


1300/5208 

NO. OF  PLY  ANGLE*  2 
ANGLE*- I 0.0 
ANGLE*  70.0 

NO. OF  LOADINGS*  2 
LOADING  UNIT  is 
•  n  riN/m 


LOADING-1 


Nl* 

4.000 

N2= 

1.000 

N6* 

0.000 

LOADING-2' 

Nl* 

l .  250 

N2= 

3.250 

N6= 

l .  299 

AFTER  5  ITERATIONS 

TOTAL  THICKNESS  * 

1 . 133E  00  cm 
90.68  PLIES 

ANGLE*- 10.0 
RATIO*  0.500 

*  PLIES*  45.34 

ANGLE*  70.0 
RATIO*  0.500 

*  PLIES*  45.34 

5’ftENGTH  RATIOS 

!*ULT I MATE  STRAiN 
AND  > I *SAF£ 

LOADING-1 
FLY  RAT IC 

:.288 

72.2  \.222 

_  I  A_i  I N  S  -  2 
3k  f  .RATIO 

-!2.0  1.068 

**  T»  ^ 

c  •  L,  .  .  -  • 


T300/5208 

NO. OF  PLY  ANGLE*  4 
ANGLE*  0.0 
ANGLE*  45.0  • 
ANGLE*-45. 0 
ANGLE*  90.0 

NO. OF  LOADINGS*  4 
LOADING  UNIT  is 
i n  MN/m 


LOADING-1 


Nl* 

1.000 

N2= 

1.000 

N6* 

1.000 

LOADING-2 

Nl* 

2.000 

N2* 

1.000 

N6= 

0.000 

LOADING-3 

Nl* 

2.000 

N2* 

-1.000 

N6* 

0.000 

LOADING- 4 

Nl* 

0.000 

N2= 

0.000 

N6* 

-1.500 

AFTER  4  ITERATIONS 

TOTAL  THICKNESS  = 
7. 727E-01  cm 
61.82  PLIES 

ANGLE*  0.0 
RATIO*  0.301 
ft  PLIES*  18.62 

ANGLE*  45.0 
RATIO*  0.254 
4  PLIES*  15.72 

ANGLE*- 45.0 
RATIO*  0.232 
#  PLIES*  14.35 


Copy  available  to  DTIC  does  sot 
BWWit  fully  legible  lepioduodos 


ANGLE*  90.0 
RATIO*  0.212 
0  PLIES*  13.12 


STRENGTH  RATIOS 

1=ULTIMATE  STRAIN 
AND  >1*SAFE 


LOADING-1 


PLY 

RATIO 

0.0 

1.320 

45.0 

2.271 

-45.0 

1.074 

90.0 

1.384 

LOADING- 

2 

PLY 

RATIO 

0.0 

1.839 

45.0 

1.476 

-45.0 

1.513 

90.0 

1.298 

LOADING- 

3 

PLY 

RATIO 

0.0 

2.396 

45.0 

1.260 

-45.0 

1.264 

90.0 

1.000 

LOADING- 

4 

PLY 

RATIO 

0.0 

1.269 

45.0 

1.010 

-45.0 

2.453 

90.0 

1.271 
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T300A5208 

BORON/5505 

AS/3501 

NO. OF  PLY  ANGLE*  4 
ANGLE*  0.0 

ANGLE*  00.0 

ANGLE*  45.0 
ANGLE=-45. 0 

NO. OF  PLY  ANGLE*  4 
ANGLE*  0.0 

ANGLE*  90.0 

ANGLE*  45.0 
ANGLE=-45.0 

NO. OF  PLY  ANGLE*  4 
ANGLE*  0.0 
ANGLE*  90.0 
ANGLE*  45.0 
ANGLE*- 45.0 

NO. OF  LOADINGS*  2 
LOADING  UNIT  is 
i n  fIN/m 

NO. OF  LOAOINGS*  2 
LOADING  UNIT  is 
MN/m 

NO. OF  LOADINGS*  2 
LOADING  UNIT  is 
!1N/m 

LOADING-1 

Nl*  3.000 
N2=  2.000 
N6=  0.500 

LOADING-1 

NI*  3.000 
N2=  2.000 
N6=  0.500 

LOADING-1 

Nl*  3.000 
N2=  2.000 
N6=  0.500 

LOADING-2 

Nl*  1.000 
N2=  4.000 
N6=  0.000 

LOADING-2 

Nl*  1.000 
N2=  4.000 
N6*  0.000 

LOADING-2 

Nl*  1,000 
N2=  4.000 
N6=  0.000 

AFTER  3  ITERATIONS 

AFTER  5  ITERATIONS 

AFTER  7  ITERATIONS 

TOTAL  THICKNESS  * 
l. 154E  00  cm 

92.32  PLIES 

TOTAL  THICKNESS  * 

1.  194E  00  cm 
95.56  PLIES 

TOTAL  THICKNESS  = 

1 . 000E  00  cm 
80.04  PLIES 

ANGLE=  0.0 

RATI  0-  0.203 

4  PLIES*  18. ?5 

ANGLE*  0.0 

RATIO*  0.203 

4  PLIES*  19.41 

ANGLE*  0.0 

RATIO*  0.203 

4  PLIES*  16.26 

ANGLE*  90.0 

RAT  1 0=  0.335 

4  PLIES*  30.97 

ANGLE*  90.0 

RATIO*  0.337 

4  PLIES*  32.2! 

ANGLE*  90.0 

RATIO*  0.333 

4  PLIES*  26.67 

ANGLE*  45.0 

RAT 10=  0.230 

4  PLIES*  21. 30 

ANGLE*  45.0 

RATIO*  0.229 

4  PLIES*  21,97 

ANGLE*  45.0 

RATIO*  0.23! 

4  PLIES*  18.56 

ANGLE=-45. 0 

RATIO*  0.230 

4  PLIES*  21.30 

ANGLE*- 45. 0 

RATIO*  0.229 

4  PLIES*  21.97 

ANGLE* -45. 0 

RATIO*  0.231 

4  PLIES*  18.56 

STRENGTH  RATIOS 

’•“ULTIMATE  STRAIN 
AND  > l “SAFE 


LOADING- 1 


PLY 

RATIO 

0.0 

l.  426 

90.  0 

l.  004 

45.  0 

l ,  324 

-45.0 

1 . 043 

LOADING- 

2 

PLY 

RATIO 

0.0 

l.  000 

90.0 

1  .  ?69 

45.0 

1. 220 

-45.0 

l  ,  220 

STRENGTH  RAH  OS 

!=ULTIMATE  STRAIN 
AND  > 1 "SAFE 


LOADING- 1 


PLY 

RATIO 

0.0 

1.  379 

90.  0 

1. 002 

45.  0 

1.  32! 

-45.0 

1 . 046 

LOADING- 

2 

PLY 

RATIO 

0.0 

l.  000 

90.  0 

1,  654 

45.  0 

1 .  300 

-45.  0 

1  .  300 

STRENGTH  RATIOS 

I -ULTIMATE  STRAIN 
AND  >!--■  SAFE- 


LOADING-! 


PLY 

RATIO 

0.  0 

1  ,  362 

90.  0 

1.  019 

45.0 

l .  287 

-45.0 

1  ,  054 

LOADING- 

-2 

PLY 

RATIO 

0.0 

1 . 000 

90.  0 

1 .  591 

45.  0 

1  195 

-45.0 

,.195 
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APPENDIX  E 
PROGRAM  LISTING 


10:  "OPTIfl":  CLEAR 
\2‘  DIM  H<61,R<31, 
S<3),  T<6>,  W(6) 
«  X(61>  Y(3 ),  Z(6 
) 

14;  OiM  A(3,  3),  0(6 
,  9),  C<6,  9),  0(3 
,  3).  £(4,  3),  G(  3 
,  3>,  N<4,  3 > »  p ( 3 
,  3>,  0(3,  35 
16; DIM  A$(6),  C#( 1 
0,  2) 

10:  l  T=i ;  (ft -10 
20.OATA  VI",  2'-, 
3",  n4M,  ,;5':,  "6!: 
21; RESTORE  20. FOR 
l  =  1  TO  6. READ  A 
*( I ) ; NEXT  i 
30; GQSU8  1200 
40. GOSUB  1700. 
GOSUB  1 100; 
GOSUB  1000 
42. GOSUB  800 
45.  IF  F$^:,‘F”GOTO 
1300 

50; GOSUB  700 
60. ! T -  I r  f 1 
70.  IF  F$-'F'  QR  IT 
>  U1THEN  GOTO  1 
300 

75: GOTO  42 
200;G*-"p,,.NC=0 
210: FOR  P=1 TO  NP 
215: IF  H(P)=0GQTO 
280 

220.  U  -P- GOSUB  650 
225; FOR  N=1 TO  NL 
230; CN=-1 
235- FOR  K*l  TO  3. 

FOR  J -1  TO  3 
240: CN=CN+G(k,  J)*£ 
( N,  J  5  *E  ( N,  K  )  . 
NEXT  / 

245-CN-CN+S(k)*E(N 
,  K )  • NEXT  K 
250. IF  CN>0LET  G*= 
■•f-'.GOTO  295 
253; IF  CN<- . 1  GOTO 
275 

255; NC*NC+ I 
260;C»(NC,  t)-CHR$ 
p; C*(NC, 2)- 
CHR$  N 


275. NEXT  N 
280; NEXT  P 
295. RETURN 
300; NM-0,  1 1  -P; 

GOSUB  650 
305. FOR  L=1 TO  NP 
310; IF  H(L)-0GOTO 
370 

315. It =L. GOSUB  600 
320. FOR  J -1  TO  3- R( 
J  1 -0 

325; FOR  k=lTO  3 
330 .  R(  J  5-R (  J  >-Q(  J, 
k)*E(N,.k) 

335 .'NEXT  k.  NEXT  J 
340; FOR  JM TO  3.V( 
J  >-'0 

342. FOR  k - 1 T 0  3 -  V ( 
J)-Y( J )+P( J, kl 
*R(K> 

345. NEXT  k. NEXT  J 
350;  ui(  L  )=0 
352-  FOR  JM  TO  3. 

FOR  kM  TO  3 
354.U(L)*W(L>+G(J. 
k>*(Y(J  i*E(N,  k 
l+£(N,  Jl*Y(k>) 
356. NEXT  k 
358.U(L)-I4(L)  +  S(  J) 
*Y( J  > • NEXT  J 
360;Nn-NHt-w(u*w(L 
) 

370. NEXT  L 
375.Nn=T(Nn> 

380. fqr  L- 1  TO  NP 
385.  uj(L>-W(L>'Nn. 

NEXT  L 
390; RETURN 
400:  FOR  1=1  TO  3- 
FOR  JM TO  3 
*05. E(  I,  J>=A(  I,  ji* 
0<  l,  J  1*S 

410. NEXT  J  NEXT  l 
415.  T=£(  i,  I  ),  U-E (  1 
,  2  >>  U-E (  1 ,  3>,  uJ 
=E(2,  2),  X  =  E(2, 
3>,  y=E<  3.  3> 
420.  QM*W*y+2*U*Xii 
U-W*U*U-T  *X*X- 

v*u*U 

425,  F  =  <t4*v-x*x>'D 
Qa(T*V-IJ*U)''D 
Ga(U*X-U*V)/0 


430.R=(T*w-u*U)/'D; 
H=.-(U*X-W*U>/D. 
Q=(U*U-T*X)/D 
435.  P(  1,  11=F,  P(  1.2 
)«G,  P<1.  3>»H,  P 
(2,  1  >*G»  P(2,  2) 
-0 

440.  p(2,  31 -Q,  P(  3.  1 
>-H,  P(3.  2)=Q,  p 
(  3,  31 -R 

450. FOR  1=1  TO  NL- 
FOR  J  =  1  TO  3 
455.  E(  I,  J  l -0 
460. for  k  —  1 T Q  3 
470.  E  (  I,  J  ) *E (  I •  J1  + 
P(  J,  k;*N(  I,  kl 
480; NEXT  k.NEXT  J 
490. NEXT  I  RETURN 
500.  FOR  1*1  TO  3- 
FOR  JMTO  3 
510;  A( I,  J  >  —0 • 0 < 1 > J 
i=0 

520. NEXT  J : NEXT  I 
530. FOR  1=1  TO  NP-  I 
1=1 

540. GOSUB  600 
550.  FOR  J--1T0  3- 
FOR  kM  to  3 
560.  A(  J,  k  >*A(  J,  k>  + 
Q(J,  k)*H< I ) . Q( 
J,  k)=0(  J,  k)+Q( 
k>*z<  i  > 

570. NEXT  k.NEXT  J. 

NEXT  I 
580. RETURN 
600. fqr  J=ir0  3 
FOR  kM  TO  3 
610. 0(  J.  kl--C(  I  I,  J* 
k) 

620. NEXT  k.NEXT  / 
640. RETURN 
550. FQR  J=ito  3 
FCR  k*]TO  3 
660. G(  J.  k)-B(  I  I,  J* 
k) 

670. NEXT  k.NEXT  J 
680.  S<  1  >  =  B<  If,  7) .  5. 
<  2  >=B<  1  f,  81  S< 
31  =B<  1  1,  5) 

690. RETURN 
700.SHME10 
732. FOR  1=1  TO  NP 
734. IF  Z< I  X >0LET 
S=-H(  I  W<  !  > 
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rr-.r'  ■J’.-LP;'  iP.’?7 


736: IF  S>0ANO  S<SM 
THEN  LET  sn=s 
737: NEXT  I 

•)1Q'  C%-"P" 

740: IF  SM>10LET  F* 
*"F": RETURN 
744:  S1=0: S2-SM: S=S 
n 

746:  IF  NO0GOTO  77 
5 

748: GOSUB  400: 

COSUB  200 
750:  IF  G*=I’F"LET  S 
2=S 

752:  IF  G*=,!P"L£T  S 
t-S 

754: IF  Sl=SnGOTO  7 
65 

755: S=<Sl+S2)/2 
756.  IF  (S2-S1X1E- 
5AND  SI “0 THEN 
LET  F*--"F" 

758:  IF  F$="F"GOTO 
786 

760:  IF  S1/CS2-S1X 
<*GOTQ  748 
762:S=SW2 
765-SR=0 

767: FOR  1  =  1  TO  NP.H 
U>=HU)+ZU>* 
S 

768:  IF  HC  I  X  IE-5 
LET  HCD=0 
770:SR=SR+HU)*HU 
> .  NEXT  l 

775- S=0:COSUB  500- 
GOSUB  «00 
778; LET  $R=T<SR) 
780. GOSUB  900 
782:  IF  CSR-SX1E-5 
let 

784;  F" OR  I  - 1  TO  NP  H 
U.)=HU>*S/SR. 
NEXT  I 

786: S=0- GOSUB  500: 

GOSUB  400 
790: GOSUB  200: 

RETURN 
800: NN=1 «  2=0 
802.  FOR  I  - .1  TO  NP.  X 
( l)-0 

804:  Z=Z*-SGN  HU). 

NEXT  I 
806. Z-l'fZ 
807. IF  NC=0GOTO  85 
0 

808. FOR  1*1  TO  NC 


7T- 


810: P=ASC  C$i I, 1 >: 
N=ASC  C*CI» 2). 
GOSUB  300 

812: FOR  J  =  i  TO  NP*  X 
( J  )=X(J  )-W.(J) 
814; NEXT  J ; NEXT  1 
820:Nf1=0:  FOR  J*lTO 
NP 

825:NM=Nf1+X<J.>*X(J 
) : NEXT  J 

830:Nn=T(NFI):  TS=0 
850:  FOR  1*1  TO  NP.'X 
( i  )=xu  )/Nri 
855: TS=TS+X< I )*Z * 
SGN  HC I ):NEXT 
l 

860:Nf1=0:FOR  I-1T0 
NP 

865: Z< I >=X(I)-TS*Z 
*SGN  HU) 

870;Nn=NI1+Z(  I  )*ZC  1 
) . NEXT  1 

ay?  P*rr‘!P'' 

874: IF  Nd< 1E-6LET 
F$=‘F,: 

876:  IF  F$*"F" THEN 
RETURN 
878.  Nfl-T  CNfl) 

880: FOR  1*1  TO  NP: Z 

<  d=z<  l  >/NH; 

NEXT  I 

885; GOSUB  500 
890: RETURN 
900. FOR  P* .1  T0  NP 
905: IF  H(P)=0GOTQ 
980 

910. I  I =P: GOSUB  650 
915. FOR  N« 1  TO  NL • 8 
=0; C=0 

920; FOR  1*1  TO  3- 
FOR  J-.1  TO  3 
925.C=C~SR*SR*GU, 
J>*E<N,  l  )*ECN, 
J) 

930. NEXT  J 
940. 8=B-SR*Sf I >*£< 
N,  I  ) .  NtX  T  I 
950: SU=C-B+T(B*S-4 
*C*C1 -  IE-6) ) )' 
'2XU-1E-6)  > 
960: IF  SU>SLET  S=S 
U 

970: NEXT  N 
980: NEXT  P: RETURN 

Odd.  Z  -  ’  'S  <  NP  > 

' Z 1 0 • r OR  I  =  >T0  NP 
Z(  I  )=Z. HC I ) 
=Z;next  1 


1020. GOSUB  500 
1030: S=0. SR= 1 
GOSUB  400 
1040, S=0. GOSUB  90 
0 

1050.  for  I -.  iro  NP 
HC I ) -H <  I  )*S 
NEXT  I. 

1060 • S=0 ■ GOSUB  50 
0 

1070. GOSUB  400. 
GOSUB  200. 
RETURN 

1100:  FOR  I  -  .!  TO  NP 
1110; C2=C0S  ( 2*T C 
l  > ) : C4=C0S  < 
4 *TU  > ) 

It 20: S2=SIN  <2*TC 
I ) ) :  S4--SIN  ( 
4*TU  )) 

II 30.  BC  I,  1  XUUC2 
*U2  +  C4iU3.  BC 
I,  -C2*U 

2+C4*U3 

1140.  BU,  2 )  =U4 -C4 
*U3-  BC  I,  9 ) -U 

5- C<**U3 

1150;  BU,  3  )=S2/'2* 
02 +S4*U3- B< I 
, 6) =S2x2*U2- 
S4*U3 

1.160.  BU,  ?)*06+C2 
*U7 :  B  <  1 ,  8  )=U 

6- C2*U7:  BU, 
5 )=S2*U7 

1 170: CC1,  1 >*U1  +  C2 
*U2+C4*U3: C( 
I,  4>--ui-C2*U 
2+C4TU3 

1180:  C Cl,  2>=U4-C4 
*U3-  CC  I,  9XU 
5-C4TU3 

1 1 90: CC 1,  3)=S2/2* 
U2+S4TU3- CCI 
,  6)=S2''2*U2- 
S4*U3 

M 95; NEXT  !  . 
RETURN 

’200. GOSUB  1600. 

LF  Cl) 

1210; GOSUB  1830 
1220: INPUT  "NO. OF 
PLY  ANGLE'S 
p"  ;  NP 

'  230: LPRINT  "NO.  Q 
f  PLY  ANGLE* 
"iNP 

' 235 • -OR  :*• “0  NP 


18 


Copy  available  to  DTIC  does  BOl 
permit  fully  legible  reproduction 


1240; PAUSE  "PLY  A 
NGL£-"+STR* 
<1): INPUT  T < 
l) 

124]; GOSUS  1810 
1242:LPRINT  "  AN 
GLE=’:;  T(I ) 
1245: NEXT  l : LF  (1 
) 

.1248:  GQSUB  1830 
1250: INPUT  "NO. OF 
INOEP.  LOAD 
IN6S=" » NL 
1252: LPRINT  "NO.Q 
F  LOADINGS=“ 

;nl 

1253; LPRINT  "  L 
OADING  UNIT 

•  $  in  !1 

N/m  "  •  LF  (1 
> 

.1255.  FOR  1=1  TO  NL 
: LPRINT  "  L 
OADING-"+A*< 
I) 

1256: FOR  J=1 TO  3 
’260: L=J: IF  J-3 
LET  L=6 

1262: X*="  N"+A 

♦  C  L  )  +  "  =  " 

1265: GOSUB  1800 
1267: LPRINT  X*. 

INPUT  NCI,  J> 

: LF  <-!> 

1268: LPRINT  NCI, J 

1269:  NCI,  J)-N(l,  J 
)*1E6 

1270:  NEXT  J • LF  (1 
> : NEXT  I . 
RETURN 
1300: ts=0 
1310: for  1  =  1  to  HP 
1320: TS=TS+H(I): 
NEXT  l 

1322.  LPRINT  AF TE 
R  " +STR$  (IT 
>+"  I TERATIQ 
NS"  •  LF  <1.) 
1325. GOSUB  ! 840 
1330. LPRINT  -TOTA 
L  THICKNESS 
=  ", T 5*180- LF 
(-1.) 

1340:  LPRINT 

cm" 

1345:  B=INT  (CTS/T 
P+.005)*100) 
'188 


1347: GOSUB  1850 
1350: LPRINT  -  " ; 8 
PLIES" 

1355: GOSUB  1500 
1360: GOSUB  1400: 

GOTO  20 
1400; LF  (1) 

1405 : LPR I  NT  "  STR 
ENGTH  RATIOS 
".LF  Cl) 

1410: LPRINT  "  1 =U 
LTIMATE  STRA 
IN  AND  >1 
=SAFE" . LF  (1 

1420: FOR  1=1  TO  NL 
1430: LPRINT  » 
LOADING-11-* 
STR*  Cl) 

1440: LPRINT  "  PL 
Y":  LF  (-1) 
1445: LPRINT  » 

RAT 

10" 

1446: GOSUB  1450 
1448: NEXT  I: 
RETURN 

1450: FOR  P=1T0  NP 
1452: IF  HCP>=0 
GOTO  1490 
1455: II =P: GOSUB  6 
50 

1460: A=0: B=0 
1465: FOR  J=ITO  3: 

FOR  K=1T0  3 
1470: A=A+GCJ,  K)*E 
Cl,  J)*ECI,  K) 
: NEXT  K 

1475:0=B+SCJ)*ECI 
,  J) : NEXT  J 
1480:  A=C-B-*-«T  CBKB+ 
4*A))/2/A 
1482: A=INT  CCA+.5 
/1E3)*1E3V1 
E3 

1483: GOSUB  1810 
1484: LPRINT  "  "J14 
CP) : LF  C-l) 
1485: GOSUB  1800 
1486: LPRINT  A 
1490: NEXT  P 
1495: LF  Cl): 
RETURN 

1500: FOR  1  =  1  TO  NP 
1505: LF  Cl) 

1510: WCI )=INT  CCT 
CI)+.05)tl0) 
/\0 

1520: A-INT  CCHCI) 
/TS+.5/1E4)* 
1E4)/1E4 


t530-  B=.INT  CCHCI) 
'TP* . 005)*10 
0V100 

1535. GOSUB  1810 
1540. LPRINT  "ANGL 
E=";ul(I  ) 

1545. GOSUB  1800 
1550: LPRINT  "RATI 

o=" ;  a 

1555. GOSUB  1850 
1560: LPRINT  "4  PL 
IES-" ;  8 
1570: NEXT  I ; 

RETURN 

1600: INPUT  HATER 
IAL  NO.  -■  ",f1 

* 

1602-  l=UAL  11**10+ 
1600: GOSUB  I 
: RESTORE  I. 
READ  Ul$,  Q,  0, 
R,  T,  Tp,  u,  U,  kl 
,  X,  V;  LPRINT 
14*.  GOTO  1697 
1610: DATA  "T300/5 
208", 181 , , 10 
3,  .28,  7.  17, 
125E-6,  1500, 
1500,  ^0. ,  2-»6 
,  68.  RETURN 
1620: DATA  "BORON/ 
5505", 204., 1 
8.5,  .23,  5.59 
.  125E-6,  1260 
.  2500,  61 .  ,  .20 
2.,  67.  . 
RETURN 

1630. DATA  AS/350 
1'  ,  138.,  8.96 
.  . 30,  7.  ! ,  125 
E-6,  1447.  144 
7,  51.  7,  206., 
93, ■ RETURN 
1640. DATA  "  ' 
RETURN 

'650.  DATA  ' 
RETURN 
1660. DATA  " 

RETURN 
1670. DATA  " 

RETURN 

1680v  DATA  "  ' 
RETURN 
1690. DATA  " 

RETURN 
1597. RETURN 


klA'.  v\  o'o  O 
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1700:  0-0*1  E9.  Q=Q* 
IE9>  t=T*1£9, 
U-U*1£6i U=U* 
IE  6*1  W=W*lE6i 
X=X*1£6,  V-Y* 
1E6 

1 705 : NN= 1/(1 -R*R* 
0/0): 01 -NN*0 
■Q2-NH*Q:Q3= 
NN*R*Q: Q4  =  T 
1710: Ul~ <  3*Q1+3*Q 
.2+2*03+4*04) 
'8- U2=<  Q1-Q2 
)/2:U3-<Ql+Q 
2-2*03-4*04) 
'8 

1720: U4=(01+Q2+6* 
Q3-4*Q4)/8:U 
5=<Q1+Q2-2*Q 
.  3+4*Q4)/8 

1730:  Q.=  1 / ( U*U ) : Q= 

1  v ( id*X ) :  Trt ./ 

( Y*Y) : G=l/U- 

1  /O:  H=.l  /U-1/ 
X 

1740:  (U*Q)/'2 

1750: GP=0*Q 1*0 1+2 
*F*Q1*Q3+Q*Q 
3*03- GQ-Q*Q3 
*Q3+2*F*Q3*0 
2+Q*Q2*Q2 
1760:  GR=0 *0.1  *Q3+F 
*(Q1*Q2+Q3*Q 
35+Q*Q3*Q2.G 
S=T*Q4*Q4 
1765. GX=G*Q1+H*Q3 
: GV=G*Q3+H*Q 

2 

1770- U1~(3XGP+3*G 
Q+2*GR+4*GS) 
'8.U2=(GP-GQ 
)/2 

t775:  U3~<GP*-GQ-2* 
GR-4*GS )/8-U 
4=<GP+GQ+6*G 
R-4*GS>'8 
1780.  U5=(GP+GQ-2* 
GR+4*GS)/8 
l  ?90- U6-CGX+GY  V2 
U7-<  GX-GV )/ 
2. RETURN 

1800. USING  '*»*»♦».*» 
4 4 ".RE  TURN 
1810:  USING  "444-.4 
"•.RETURN 


1830* USING  "44": 
RETURN 

1840: USING  "4444. 

444a". RETURN 
1850- USING  "4444. 

44". RETURN 
1890: END 
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